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We measured the coherence properties of the free-electron laser SACLA on a single shot basis at an X-ray 
energy of 8 keV. By analysing small-angle X-ray scattering speckle patterns from colloidal dispersions we 
found a degree of transverse coherence of ^, — 0.79 ± 0.09. Taking detector properties into account, we 
developed a simulation model in oder to determine the degree of coherence from intensity histograms. 
Finally we calculated a coherence time of t^— 0.1 fs and a pulse duration of 5.2 fs which corresponds with 
previous predictions. 

The rise of X-ray free electron lasers (XFEL) allows the study of physical and chemical properties that could 
not be accessed experimentally before. In particular, the ultrashort pulse duration on the order of 10 fs at 
hard x-ray wavelengths of around 1.5 A offers the possibility to access physical and chemical processes on 
molecular time and length scales by single shot x-ray diffraction patterns. 

Many of such studies rely on the coherence properties of the X-ray beam. This includes, e.g., the study of 
dynamic properties via X-ray photon correlation spectroscopy'"'', imaging single particles and biological speci- 
mens by coherent diffractive imaging'''^, cross correlation techniques to reveal a single particle's shape*" and 
orientational order in liquids and glasses', or the imaging of ultrafast magnetization processes'*. XFEL sources 
such as the Linac Coherent Light Source (LCLS) at SLAC and SACLA at SPring-8'° are based on self-amplified 
spontaneous emission (SASE)" ". The SASE process originates from random fluctuation in the electron beam. 
Therefore, the coherence properties of such an X-ray source differ from optical lasers. In particular, every single 
X-ray shot is expected to show certain random features because of the stochastic nature of the SASE process. This 
was observed recently in experiments at the LCLS and in simulations on the longitudinal coherence properties, 
where a significant shot-to-shot variation has been observed". In contrast, a SASE beam is expected to show a 
superior transverse coherence'^. 

Recently, several studies have been performed on the different EEL sources' coherence properties. In the soft x- 
ray regime, pinhole and double slit experiments have been performed at FLASH'* '"" and LCLS"" in order to study 
the coherence lengths and times. In further studies, interferometry methods have been used"''* to determine 
longitudinal and transverse coherence lengths. In contrast, coherence properties in the range of hard X-ray have 
been examined by analysing coherent diffraction patterns, so-called speckle patterns, from disordered samples on 
a shot-to-shot basis'"* ''^. A similar study has been performed at FLASH™. In this study, the degree of coherence is 
obtained from the speckle contrast. By analysing single shot diffraction patterns from colloidal particles in small 
angle X-ray scattering (SAXS) geometry, the degree of transverse coherence was found to be almost 1"'". 
Furthermore, the longitudinal coherence was determined by wide angle scattering (WAXS) experiments on 
amorphous metals and it also appeared to be very sensitive on the SASE fluctuations^ '". 

Here we present a study of the coherence properties of the free-electron laser SACLA in the hard X-ray regime. 
In contrast to the LCLS, three bunch compressors and short-period in-vacuum undulators are used at SACLA 
that result in a compact design compared to LCLS'°. Furthermore, initial experiments reported pulse durations 
well below 10 fs^'. The question arises if and how this design combined with the short pulse duration have an 
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impact on the coherence properties of the single x-ray pulses. We 
calculated the degree of coherence for single shots from speckle 
patterns of colloidal particles at a hard X-ray energy of 8 keV. The 
number of transverse modes was determined by the need of correct- 
ing for the energy bandwidth of the x-ray pulses. The results are 
supported by simulations of the impact of the charge sharing process 
of CCD detectors on intensity histograms of speckle patterns. Finally, 
we estimate the pulse duration from the intensity fluctuations of the 
beam and compare it to the expected values in the order of 5 fs. 

Results 

A typical speckle pattern from an amorphous sample together with 
the azimuthally averaged intensity I{q) is shown in Fig. 1 . The experi- 
mental details are given in the Methods section. The size of individual 
speckles depends on set-up parameters only and is given by 5 = M/l, 
with wavelength 1, sample-detector distance d and beam size at the 
sample /. The experimental speckle size can be determined by cal- 
culating the spatial intensity correlation function g(R) = {I{r + 
R)I('')}/{I(^)y-> that is connected to the coherence properties of the 
beam via the mutual coherence function y{R, z) by ^(i^) = |y(i?, z = 
0)|^ + 1 at the spatial distance -R and time difference z^^. yiR, z) is a 
complex function that contains all information of the coherence 
properties. The degree of coherence for a single speckle pattern is 

given by c= |y(i^ = 0,T = 0)| = 1 y/ Vm, with M denoting the mode 

number of the beam"'^'. By modelling g(R) with a Lorentzian func- 
tion for both the vertical and horizontal directions, we obtain a beam 
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Figure 1 | (a) Single shot speckle pattern and corresponding I{q) averaged 
over 500 shots (inset), (b) Typical intensity autocorrelation g{K) and 
Lorentzian fit to determine the speckle size. Due to the finite size of the 
extracted region of interest no full decorrelation (g(R) --1=0) can be 
observed. 



size of(1.9±0.1)X(1.6±0.1) pm^ (h X v) which corresponds with 
the directly measured beam size of 1.8 X 1.5 |im^ excellently. 

By analysing the speckle pattern the most direct quantity that can 
be accessed is the speckle contrast or visibility [i^. It is given by the 
normalized standard deviation of the intensity yielding^^"" 



(I) 
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Due to the equality j)'(i^ = 0,T = 0)j — 1^ vM, the speckle contrast 

directly reflects the degree of coherence of the x-ray beam. A fully 
coherent beam would have a mode number of M = 1 and thus a 
contrast of /i^ = 1. The experimentally measured contrast is affected 
by both the longitudinal and transverse coherence of the x-ray beam 
yielding /i^ = PtPiiq), with the transverse coherence fi^ and a q- 
dependent correction Piiq) that depends on the energy bandwidth 
and the scattering geometry"'^^'^''. 

In Fig. 2 (a) the contrast /i^ is shown for 2000 single shots atq = 
0.041 ± 0.001 nm"'. The average scattering intensity in this q-range 
was 35 photons/pixel. We find an averaged contrast of (fSs) = 0.70 ± 
0.08 which can be taken from the histogram of the contrast values as 
shown in Fig. 2 (b). This corresponds to a mode number of M = 2.0 
± 0.5. The single shot contrast varies from slightly above 0.5 up to 
some fully coherent shots with contrast values of 1 following an 
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Figure 2 | (a) Single shot speckle contrast from 2000 consecutive shots 
calculated atq = 0.041 nm"'. (b) Histogram of /Jj. The solid line represents 
a lognormal distribution using a mean value of 0.7 and a width of 0. 105. (c) 
Speckle contrast as function of q averaged over 2000 consecutive single 
shots. The error bars represent the standard deviation of the single shot 
values from the averaged value. 
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asymmetric distribution. This is demonstrated by the soUd Hne in 
Fig. 2 (b) reflecting a lognormal distribution. In addition, we did not 
find any significant variation as a function of q, see Fig. 2 (c). Taking 
into account the natural bandwidth of the FEL radiation of Aa/A = 5 
X 10"', the correction accounting for the effect of the longitudinal 
coherence can be calculated to""'^ 



dx(L — x) 



dxw-y) 
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[exp( - 2| Ax + B^t^l ) + exp( - 2| Ax + Byl )] 
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studied here, we find a correction factor of /i/ = 0.89, i.e. the trans- 
verse coherence is found to be Pt = t^JPi — 0-79 ± 0.09. This corre- 
sponds to a transverse mode number of = 1.6 ± 0.4. 

Typically, a certain minimum intensity is necessary for contrast 
calculation using Eq. 1. At low scattering intensities the contrast is 
usually obtained by describing the intensity distribution by intensity 
probability distribution functions^'^"'^'. This method was shown to 
sufficiently describe the contrast of low-intensity speckle patterns in 
particular at large q ^ and by using split-and-delay devices^'. Due to 
charge sharing of neighbouring pixels in the CCD detector and read 
out noise a continuous histogram is measured with peaks at intensity 
values corresponding to multiples of single photon events rather than 
discrete photon events. Typically, experimental intensity values are 
discretized by rounding up or down to integer values. This may intro- 
duce some artifacts. As the ratio of single, double, triple, . . . photon 
events is crucial for obtaining the mode number'"*, the interpretation 
can have a strong impact on the resulting contrast. At low intensities, 
such problems can be avoided by using so-called droplet algorithms 
to identify the location of the photons impacting the pixels of the 
detector*. However, those algorithms are limited to very low intensities, 
in particular work best for non-overlapping single photon events^". 

Therefore, here we use the experimental, non-discretized intensity 
histograms and model those by simulations of intensity patterns 
following a negative binomal distribution and take into account 
charge sharing effects of the CCD detector^^ ™ and detector noise, 
see Methods for details. We used M = 2, which corresponds to the 
experimental value discussed above, as input parameter for the mode 
numbers. The results are shown as red line in Fig. 3. For comparison, 
the black circles represent the discrete values before blurring. In 
particular, at low intensities the results correspond to the experi- 
mental data. At higher intensities (see Fig. 3 (b)), the model describes 
the general shape of the experimentally obtained histogram very well, 
however, small differences can be found for I i 0.5 and also at higher 
intensities. We think this is due to the simplicity of our model and 
due to a better statistic in the simulation model. Because of blurring, 
real scattering intensities can appear as lower intensities. For 
example, a single event of intensity 7o hitting the CCD sensor in 
between two pixels, will appear as two independent events of intens- 
ity IJl. These can not be distinguished from other events of intensity 
I0/2 in the same pixels. As a result, experimentally obtained mode 
numbers from intensity histograms may be overestimated. 

Finally we want to discuss the impact of the shot- to-shot intensity 
distribution. Thanks to the energy bandwidth of 40 eV, we estimate 
the coherence time = ^//c = 0.1 fs, with the speed of light c. The 
coherence time is proportional to the average x-ray pulse duration 
(tp) via (tp) = MiZ^", with the average number of temporal modes M;. 
We determined M; from the shot-to-shot intensity distribution via a 
Gamma-distribution""''' 
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Figure 3 | Intensity histograms at q = 0.071 ± 0.003 nm"' (a) and q = 
0.053 ± 0.003 nm"' (b) averaged over 600 consecutive shots and 
simulation results taking blurring and noise into account. 

This equation is valid at least for the linear FEL regime"*^. The pulse 
intensity was measured by a diode located in front of the optical 
components of the beamline for more than 10000 single shots. The 
shot-to-shot variation of the beam intensity is displayed in Fig. 4 (a). 
By fitting its distribution with Eq. 3 we obtain M; = 52.2 ± 2.0, see 
Fig. 4 (b). This results in an average x-ray pulse length of tp = 5.2 ± 
0.2 fs which perfectly corresponds to other recent results from 
SACLA"'. 



Discussion 

In conclusion, we present a study on the coherence properties of the 
hard x-ray free-electron laser SACLA based on the analysis of single 
shot speckle patterns from colloidal particles. The size of the speck- 
les in horizontal and vertical directions are in perfect agreement 
with expected values. We found a degree of transverse coherence 
of 0.79 ± 0.09 corresponding to a mode number of M = 1.6. The 
degree of transverse coherence is affected by set-up and machine 
parameters such as bunch compression schemes, saturation length, 
and photon wavelength resulting in observed values below 1. This is 
in agreement with expectations from FEL theory"'^' and previous 
experiments""'"'™. 

By using a simple simulation model, we were able to model experi- 
mental intensity histograms. Because of the model's simplicity, the 
illustration of the values in the histogram is vague, i.e. they are 
blurred. This blurring's origin is related to the charge sharing and 
the noise level which is present in CCD detectors used in modern x- 
ray scattering experiments. Such simulation methods wUl help to 
determine contrasts more efficiently, in particular for weak scattering 
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Figure 4 | (a) Shot to shot variation of the beam intensity, (b) Intensity 
distribution. The solid line is a fit of Eq. 3 yielding M; = 52.2 ± 2, the mean 
intensity was (J) = 300 nJ. 

intensities where time-consuming droplet algorithms are used to 
overcome the effects of charge sharing'. However, we want to note 
that the simple simulation model discussed here does not give a full 
description of such influences. Further knowledge of detector prop- 
erties is needed for a full imderstanding of our observations. 

Methods 

Sample. We synthesized the sihca colloidal particles by using a modified Stober 
method^^-^*. The particles were dispersed in water at a volume fraction of about 5%. 
They had a mean radius of 41 nm with a polydispersity of 10% which was determined 
by a form factor analysis {see following section). After the preparation the dispersion 
was filled into a quartz capOlary {diameter of about 0.7 mm) which has been sealed by 
melting the glass. 

Experiment. The experiment was performed at SACLA's EH3 of BL3 at a photon 
energy of 8 keV. With a bandwidth of 5 X 10"^ at a repetition rate of 20 Hz^^, the 
pulse energy was {300 ± 40) The beam was focalized with KB mirrors to a size of 
1.8 X 1.5 ]im^ (h X v) at the sample position. The beam size was measured in the 
vertical and horizontal directions by knife edge scans. The sample capillaries were 
placed into a holder that was mounted in the MAXIC sample chamber of BL3^^. To 
resolve the structure factor of the samples, the detector {dual multiport charge 
coupled device {MPCCD}^^, i.e. two sensors with 1024 X 512 pixels each with 
individual pixel size of 50 |im X 50 |im) was placed at a distance of 3 m behind the 
sample. The accessible g-range was limited by the beam stop and the size of the 
detector to 0.02 nm~^ < q < 0.34 nm"'. For each sample at least 10000 single shot 
speckle patterns were taken. These were averaged and azimuthally integrated to 
obtain I{q) as shown in Fig. 1. We measured several 1000 speckle patterns from a 
diluted dispersion (volume fraction £0.1 %) in order to extract the particle form 
factor. This was fitted by a form factor of a sphere convoluted with a size distribution 
to account for the polydispersity^^. At the volume fraction of 5% used for the 
coherence analysis a structure factor reflectmg particle- particle correlations was 
visible^^ see Fig. 1. 

Charge sharing simulations. Depending on the detector used, the photons 
impinging on the CDD sensor generate a charge cloud of different sizes, typically in 
the order of 10 |j.m^^. This may result in shortcomings, e.g. an intensity of 1.4 can 
either be achieved by one photon with negligible charge sharing and less charge 
sharing contribution from neighbouring pixels, or e.g. by two photons (or more) with 
non- negligible charge sharing contribution. The impact of charge sharing of the CCD 
pixels and noise was estimated by computer simulations. Therefore, we used several 
500 X 500 pixel grids with values corresponding a negative binomal distribution at an 
average intensity of 0.36 and 1.54 photons per pixel {identical to the experimental 
data), respectively. The negative binomal distribution function p[,{I) is given by^^ 



where / is the intensity in number of photons, (I) its average, and T the Gamma- 
function. The initial mode number was set to M ~ 2, matching the result discussed 
above. Each pixel pattern was blurred by a random Gaussian distribution to model the 
charge sharing of pixels given by 



(x-Xc) [y-ycY 



2a^ 



(5) 



Here, x, y denote the pixel coordinate and the random center of the Gaussian. 
After that, Gaussian noise was added to the pattern. The amplitude of the charge 
sharing parameter a and the noise were adapted to match the experimentally 
obtained histograms. Therefore, we calculated an averaged histogram from 600 single 
shots showing low intensities at g — 0.071 ± 0.003 nm~\seeFig. 3 {a). The intensities 
in histograms are scaled for the ADU {Analog- Digital Unit) value for 1 photon of /{I 
photon) — 124.7 adu. By comparing the simulated histogram with the averaged 
experimental one we find ~ 0. 1 5 for the strength of charge sharing and a width of 
Gaussian noise of o"noise ~ 0.1 photons. The effect of charge sharing is only slightly 
above the expected value for the MPCCD detector of (Tmpccd S 0.10 pixels^^. 
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